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We studied the drag resistivity between dilute two-dimensional hole systems, near the apparent 
metal-insulator transition. We find the deviations from the dependence of the drag to be inde- 
pendent of layer spacing and correlated with the metalliclike behavior in the single layer resistivity, 
suggesting they both arise from the same origin. In addition, layer spacing dependence measure- 
ments suggest that while the screening properties of the system remain relatively independent of 
temperature, they weaken significantly as the carrier density is reduced. Finally, we demonstrate 
that the drag itself significantly enhances the metallic T dependence in the single layer resistivity. 

PACS numbers: 73.40.-c,71.30.+h,73.40.Kp,73.21.Ac 



I. INTRODUCTION 

Two-dimensional (2D) electron transport in semicon- 
ductor heterostructures has provided a rich venue for the 
study of electron interaction physics, exhibiting such ex- 
otic states as the fractional quantum Hall liquid-'-. Here 
strong interactions, induced by a large perpendicular 
magnetic field, stabilize this non-Fermi liquid state. Re- 
cently, much attention in the field has focused upon di- 
lute 2D systems, which are characterized by large ratios 
of carrier interaction energy to kinetic energy {ts > 10). 
These systems exhibit an anomalous metalliclike behav- 
ior and an apparent metal-insulator transition^, contra- 
dictory to the scaling theory of localization-. To date, 
the origin of the metallic behavior is unclear, with several 
fundamental questions regarding this regime still unan- 
swered. Among the most important of these are the na- 
ture of the many-body correlations and screening proper- 
ties in this dilute limit. To gain insight into both of these 
issues, we have performed frictional drag measurements. 

Drag measurement^ allows one to directly study 
carrier-carrier interactions. These experiments are per- 
formed, on double layer systems, by driving a current 
(Id) in one layer, and measuring the potential (Vd), 
which arises in the other layer due to momentum trans- 
fer. The drag resistivity (pn), given by Vd/ Id, is directly 
proportional to the interlayer carrier-carrier scattering 
rate. Furthermore, the layer spacing (d) dependence of 
the drag provides a powerful probe of both the interlayer 
correlations and the screening properties, which exist in 
the system. Any change in either of these properties 
will manifest itself in the layer spacing dependence of the 
drag. Recently, the drag was measured between low den- 
sity hole systems^, with approximately ranging from 
10 to 20^. In this regime, po showed a 2 to 3 orders 
of magnitude enhancement over the theory for weakly 
interacting systems^, and the corresponding low density 
electron results^. In addition, an anomalous temperature 
(T) dependence, which could not be explained in light of 
previously studied drag processes, was observed, with po 
exhibiting a greater than dependence at low tempera- 



tures. Upon further increase of T, a crossover to a weaker 
than dependence was observed. It was proposed that 
these deviations arose from many-body correlations in 
such a strongly interacting regime.^'^. 

In this article, we study the density, temperature, and 
layer separation dependence of the drag between strongly 
interacting 2D hole systems, in the vicinity of the appar- 
ent metal-insulator transition. Using this data, we show 
that the deviations from the dependence of the drag 
are not a result of a phonon mediated drag processiSiii, 
but rather arise from a novel mechanism related to the 
large Vg value of the system. We find these deviations 
to be independent of layer spacing, implying that they 
arise from an intralayer correlation effect. Furthermore, 
we find that the deviations to the dependence of 
the drag are correlated with the metalliclike T depen- 
dence in the single layer resistivity, suggesting that both 
anomalies have the same origin. Our layer spacing de- 
pendence data imply that while the screening properties 
in this regime remain relatively independent of temper- 
ature, they weaken significantly as the carrier density is 
reduced towards the metal-insulator transition. Finally, 
we demonstrate that the Coulomb drag effect itself can 
significantly enhance the metallic T dependence in the 
single layer resistivity. 

II. EXPERIMENTAL DETAILS 

Four different samples were used in this study. Each 
sample contains a double quantum well structure, con- 
sisting of two Si doped p-type GaAs quantum wells sep- 
arated by a pure AlAs barrier--, which was grown by 
molecular beam epitaxy on a (311) A GaAs substrate. 
Sample A, which was used in an earlier studjt^, has an 
average grown layer density and center to center layer 
separation of 2.5 x lO^*' cm~^ and 300 A, respectively. 
Samples B and C were similar to sample A, with the 
exception of having different center to center layer sep- 
arations of 225 and 450 A, respectively. Sample D was 
a higher density sample, having an average grown layer 
density and center to center layer separation of 7.0 x 10^° 
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FIG. 1: Pd/T^ vs T at p„, = 2.5 x 10^" cm-^ for diflFerent 
d. For clarity, data from the d = 450 A sample has been 
multiplied by a factor of 6.5. The dashed line marks the peak 
position. Open circles measured in the dilution refrigerator. 



FIG. 2: pD vs pm on log-log scale, at T = 300 mK. Data from 
Sample A (d = 300 A) and D (d = 275 A) are shown by soUd 
and open circles, respectively. The solid line is a fit with slope 
-5. Inset (a): pn vs pi on log-log scale, with p2 — 2.15 x 10^° 
cm~^, at T = 80 and 300 mK. Solid lines are fits with slope 
close to -2.5. Dashed line indicates matched density. Inset 
(b): PD vs Vbias- Pm = 2.1 X 10^° cm"^ and T = 230 mK. 



cm"^ and 275 A, respectively. Explicit details of the pa- 
rameters for each of these four samples are listed in Table 
I. The samples were processed allowing independent con- 
tact to each of the two layers, using a selective depletion 
schemei^. In addition, both layer densities are indepen- 
dently tunable using evaporated metallic gates. 

The data presented in this paper were obtained using 
^He and dilution refrigerators. The densities in each layer 
were determined by independently measuring Shubnikov- 
de Haas oscillations. Drive currents between 50 pA to 10 
nA were passed, in the [233] direction, through one of 
the layers, while the drag signal was measured in the 
other layer, using standard lock-in techniques at 4 Hz. 
To ensure that no spurious sources were contributing to 
our signal, all the standard consistency checks associated 
with the drag technique were performed^. 



TABLE I: Sample parameters. Mobilities quoted at 300 mK. 



Sample 



A B C D 



Average layer density [x 10^" cm"^] 2.5 2.5 2.5 7.0 

Top layer mobility [xlO"^ cmVVs] 1.5 1.9 2.9 5.7 

Bottom layer mobility [x 10^ cm VVs] 1.5 1.3 1.7 7.7 

Quantum well thickness [A] 150 150 150 175 

Barrier thickness [A] 150 75 300 100 

Center to center layer separation [A] 300 225 450 275 



III. RESULTS AND DISCUSSION 

We begin our presentation of the data by first looking 
at the temperature dependence of the drag in the dilute 
regime at different layer spacings. This is presented in 
Fig. n where we plot po/T^ vs T at matched densities 
(pm) of 2.5 x 10^° cm-^ for d = 225, 300, and 450 A. 
For clarity, the data from the d = 450 A sample has been 
multiplied by a factor of 6.5. The error bar shown for 
the d — 225 A data is due to a small interlayer leakage 
contribution. Note that both the d — 225 and 450 A 
samples exhibit the same qualitative behavior, showing a 
peak in pd/T^ vs T, as was reported earlier— at d = 300 
A. In addition, the peak position in po/T^ vs T, marked 
by the dashed line, appears to be independent of layer 
separation. Although not shown here, we also observed 
this at lower densities. 

Before discussing this data, we would like to show that 
these deviations from the dependence are not a re- 
sult of a phonon mediated drag processiSiii. It has been 
well established that a contribution from 2fci? phonon 
exchange to the drag will produce deviations from the 
expected dependence, which are qualitatively similar 
to those shown in Fig 1. However, the fractional devi- 
ation from arising from a phonon mediated process 
should decrease considerably as d is reduced. This stems 
from the fact that as d is reduced, the direct Coulomb 
component of the drag increases significantly, whereas 
the phonon mediated component increases with a much 
weaker d dependence^^. The data in Fig 1. show that 
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the fractional deviation from is roughly independent 
of d. Here, at all three layer spacings, po/T^ changes 
roughly 25 to 30 %, from 0.5 to 2.0 K. 

Furthermore, to conclusively rule out the phonon medi- 
ated process, we present the relative density dependence 
of the drag at d = 300 A, in Inset (a) of Fig. [3 Here 
we have measured po as a function of one layer den- 
sity (pi), while keeping the other layer density fixed at 
P2 = 2.15 X 10^" cm^^. It has been demonstrated that 
the phonon mediated component of the drag is strongly 
suppressed as the layer densities are mismatchediSiii. If 
2kF phonon exchange contributes to the drag in our sam- 
ples, one would expect to see some signature of a local 
maximum at matched density in these density ratio mea- 
surements. However, as shown in the inset, where pjj is 
plotted on log-log scale against pi at both T = 80 and 
300 mK, the drag exhibits no signature of a local maxi- 
mum at matched density, and is described by a linear fit 
with slope close to -2.5, implying pn oc (pip2)~^^^ at low 
temperatures^^. The deviation from the fit at low den- 
sity in the 300 mK trace is a consequence of the density 
dependent crossover in the T dependence of prA- 

To provide yet further evidence against the presence of 
2k f phonon exchange in our samples, we have performed 
density imbalance measurement^. These measurements 
are performed by measuring pjj as a function of inter- 
layer bias (Vbias)- Interlayer bias transfers carriers be- 
tween the layers, while keeping the total density in both 
layers fixed. Such a measurement provides a very sensi- 
tive probe as to the presence of phonon mediatedi2*ii (or 
direcli^) 2kF scattering. A series expansion of the power 
law (piP2)~^^^ dependence deduced from our density ra- 
tio measurements implies that po will exhibit a quadratic 
increase with interlayer bias. If po contains any signif- 
icant contribution from 2kp phonon exchange, this will 
clearly lead to the drag decreasing as a function of Vbias ■ 
In Inset (b) of Fig 2, we plot p^) vs Vbias at T = 230 
mK. Here the layer densities are matched at 2.1 x 10^° 
cm~^ at zero bias. It is clear from the figure that pjj 
exhibits a quadratic increase with Vbias, consistent with 
our density ratio results. Here T/Tp = 0.07, so ther- 
mal broadening of the Fermi surface is negligible. This 
data conclusively shows that 2kF phonon exchange^°-^^ 
(or direct 2k p scattering processes^) do not yield any 
significant contribution to the drag at c? = 300 A. 

In the main plot of Fig 2, we investigate the matched 
density dependence of po- Here vs Pm is plotted 
on log- log scale for densities ranging from pm = 1.0 to 
2.5 X 10^° cm~2, taken from Sample A (d = 300 A), and 
at Pm = 7.0 X 10^°, taken from Sample D (d = 275 A). We 
find that this data is very well described by a linear fit 
with slope of -5, implying that for this density range, cor- 
responding to rg = 5.1 to 13.5 (using m* — 0.17me), po 
follows a dependence at low temperature, consistent 
with our density ratio results. In the Coulomb drag the- 
ory for weakly interacting systems^, which is successful 
in explaining the results of drag experiments performed 
on high density electron systems^ with ~ 1, p^i is ex- 
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FIG. 3: p(T)/p(To) vs T at p = 2.5x10^°l:„. 
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pected to scale as p:^ . The stronger density dependence 
we have observed here is consistent with earlier reports^ 
that the drag between low density holes exhibits a sig- 
nificant enhancement over this simple theory, with the 
discrepancy increasing as the density is reduced. In this 
context, what is of particular interest is a recent experi- 
ment, on low density electron double layer systems with 
d = 280 A, which studied the drag in a regime of inter- 
mediate interaction strength, with ranging from about 
2 to 4.3-. This study also found that the drag exhibited 
a power law dependence on the matched electron density 
(n), and observed that po oc n^'^. This result suggests 
that there is a crossover in the behavior of the drag from 
a to a dependence as the value of the system 
is increased. This provides yet further evidence that the 
large enhancement of the drag and the deviations from 
the dependence found in the dilute regime arise from 
correlation effects. 

It is important to note that the deviations we have 
observed could arise from two different types of correla- 
tion effects: either interlayer or intralayer correlations. 
The T dependence data for different d shown in Fig 1 
provides very useful information regarding differentiat- 
ing between these two effects. If the deviations from the 
dependence arise from an interlayer correlation effect, 
we should expect them to exhibit a significant change as 
the layer spacing is varied. However, this is not what is 
observed. Both the fractional deviation from and the 
peak position are found to be independent of d, imply- 
ing that interlayer correlation effects are not playing a 
role here, and that the deviations arise from something 
intrinsic within each of the single layers themselves. 

At this point, it is appropriate to compare the devia- 
tions from the dependence of the drag to the metal- 
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liclike behavior observed in the individual layers of the 
samples used in this study. Many believe that this strong 
metalliclike temperature dependence arises from correla- 
tion effects, which exist in such a large regime. By 
independently contacting one of the layers in our double 
layer system, we can measure the temperature depen- 
dence of the single layer resistivity. However, it is im- 
portant to realize that the resistivity obtained from such 
a measurement is significantly enhanced due to the drag 
effect. Due to its large magnitude, the drag effect will 
make a significant contribution to the single layer resis- 
tivity. Since the drag has a relatively strong temperature 
dependence, this will significantly affect the strength of 
the metallic behavior obtained in a single layer resistiv- 
ity measurement. This is demonstrated in Fig. 3, where 
we have plotted the normalized T dependence of the sin- 
gle layer resistivity, dX p — 2.5 xfO^" cm^^, for different 
d. For reference, the single layer resitivities at T = 400 
mK are 1.48, 1.83, and 1.85 kVl/U for the d = 225, 300, 
and 450 A samples, respectively. While all three of the 
samples here have intrinsically similar properties (resis- 
tivity and mobility) , it is clear that the metallic behavior 
becomes much stronger as d is reduced, as should be ex- 
pected from the increasing drag contribution. Therefore, 
it is appropriate to subtract the drag resistivity from the 
measured single layer resistivity to obtain the intrinsic 
single layer resistivity- that is the resistivity which would 
be measured if the layers could be separated infinitely 
far apart. This is shown in the inset, where we sub- 
tract pd from p and plot the normalized T dependence 
of p - Pd- Note here that the metallic behavior in all 
three curves shows roughly the same strength, confirm- 
ing that the change in the metallic behavior of the single 
layer resistivity, shown in the main plot, arises from the 
drag contribution to the resistivity. The most striking 
feature here is that the temperature at which the metal- 
lic behavior is suppressed vcl p - p^ is at roughly the 
same temperature at which the peak in p^jT^ vs T is 
observed, as shown by the dashed line in the inset. In ad- 
dition, in these three samples, the metallic temperature 
dependence produces a roughly 20 to 25 % change in p - 
Py). This is extremely close to the 25 to 30 % deviation 
from T^, which p^ exhibits over the same temperature 
range. These similarities seem to suggest that the devia- 
tions from the dependence of the drag are correlated 
with the anomalous metallic behavior found in the single 
layer, and implies that they most likely have the same 
origin. 

The fact that both the single layer resistivity and the 
drag exhibit similar anomalies is quite striking, since 
these are two extremely different transport properties. 
Although there have been numerous models attempting 
to explain the origin of the metallic behavior, here we 
choose to analyze our data in light of one of the most 
prominent: temperature dependent screeningi^. In this 
model, the metallic behavior arises from temperature de- 
pendent changes in the static screening of ionized impu- 
rity potentials, which become important at large r^. As 
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FIG. 4; Pd vs d, on log-log scale, for -pra = 1.5, 2.0, and 
2.5 X 10^° cm"^ at T = 1.0 K. Solid lines are the best linear 
fits of each data set. 
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FIG. 5: a vs T for = 1.5, 2.0, and 2.5 x 10^° cm'^ a. 
deduced from fitting pn vs d to a d~° fit. 



the temperature is increased from low T, the screening 
weakens significantly, leading to the metalliclike increase 
in the resistivity. It is possible to envision that such a 
screening change could give rise to the enhancement to 
the dependence of found at low T . However, here 
we are concerned with the dynamic screening properties 
of the system as it screens the interlayer Coulomb poten- 
tial. 

If the deviations from the dependence of the drag 
arise from a temperature dependent screening effect then 
this will surely manifest itself in the layer spacing depen- 
dence of the drag. The strength of the d dependence of 
Pd is proportional to the screening strength of the in- 
terlayer Coulomb potential. If the screening is strong, 
the interlayer potential will drop rapidly with d and pu 
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will show a strong decrease as d is increased. On the 
other hand, if screening is weak the interlayer potential 
will drop off much slower with d, and po will exhibit a 
much weaker decrease as d is increased. By investigating 
the strength of the d dependence at different tempera- 
tures, we can determine whether the dynamic screening 
properties of the 2D system in this regime exhibit any 
temperature dependent behavior. Using data from Sam- 
ples A, B, and C, we have studied as a function of d at 
various temperatures, for Pm = 2.5, 2.0, and 1.5 x 10^° 
cm~^. A typical data set, obtained at T = 1.0 K and 
plotted on log-log scale, is shown in Fig 4. As expected, 
we observe that po shows a strong increase with decreas- 
ing layer spacing. The solid lines are linear fits of each 
data set, and their slopes correspond to the exponent, 
a, where po cx: d~". Using additional data at different 
temperatures, which are not shown, we can determine 
the strength of the d dependence at each density, for 
different T. We present this data in Fig 5, where we 
plot a as a function of temperature for different Pm- We 
mention that the a found for pm = 1-5 x 10^° cm~^ is 
deduced only using data at d = 225 and 300 A. Before 
discussing this data, we would first like to mention that 
due to the narrow range of data we are not making any 
quantitative claims on the values of these exponents. On 
the other hand, a comparison of the relative change of 
a with temperature and density is perfectly valid and 
independent of any of the systematic error in the d val- 
ues of our samples. The first point we make is that at 
all three densities, the exponent shows a weak increase 
with T, implying that the dynamic screening properties 
in this regime are relatively independent of temperature. 
If anything, the screening strengthens slightly as T is in- 
creased, the opposite of what should be expected from 
the T dependent screening model. Another point which 
is clear from this data is that the strength of the d de- 
pendence weakens significantly as the carrier density is 
lowered. This implies that the dynamic screening proper- 
ties of the 2D system weaken dramatically as the carrier 
density is reduced towards the metal-insulator transition. 

Finally, we would like to refer back to the single layer 
transport data presented in Fig 3, which clearly shows 
that in the dilute regime, the drag constitutes a signifi- 
cant fraction of the resistivity. In contrast, at small 
the drag contribution to the resistivity is roughly 0.1 % at 



these temperatures. This demonstrates that in the dilute 
regime an electron-electron (e-e) scattering process can 
yield a substantial T dependent correction to the resistiv- 
ity, and suggests the possibility that the metallic behavior 
seen in so many dilute 2D systems might arise from an 
e-e scattering process. For example, a theory on the in- 
teraction corrections^^ in this regime has been shown to 
qualitatively explain the metalliclike temperature depen- 
dence measured in experiments^^. In this theory the dis- 
order breaks the translational invariance of the system, 
and allows e-e scattering to make a finite contribution 
to the resistivity. Our data suggest that the magnitude 
of this interaction correction could possibly constitute a 
significant fraction of the resistivity. Another possibility 
is that the electron system forms two channels, with only 
one carrying the current, which has been suggested by 
local compressibility measurementsiS,. Here interchannel 
e-e scattering contributes to the resistance of the current 
carrying channel and could possibly create the metallic 
T dependence. This is analogous to the double layer sys- 
tem, in which only one layer carries current and the other 
layer causes dissipation via interlayer e-e scattering. 



IV. CONCLUSIONS 

In conclusion, we have shown the deviations from the 
dependence of the drag to be independent of layer 
spacing and correlated with the metallic T dependence in 
the single layer resistivity, suggesting that both anoma- 
lies have the same origin. Our studies of the strength of 
the layer spacing dependence of po imply that the dy- 
namic screening properties in this regime are relatively 
temperature independent. However, they weaken signif- 
icantly as the carrier density is reduced. Finally, in this 
dilute regime, we have demonstrated that the drag effect 
itself can significantly enhance the metallic behavior in 
the single layer resistivity. 



V. ACKNOWLEDGMENTS 

This research was funded by the NSF, a DURINT 
grant from the ONR, and the DOE. 



^ D.C. Tsui, H.L. Stormer, and A.C. Gossard, Phys. Rev. 

Lett. 48, 1559 (1982). 
^ See for a review E. Abrahams, S.V. Kravchenko, and M.P. 

Sarachik, Rev. Mod. Phys. 73, 251 (2001). 
^ E. Abrahams, P.W. Anderson, D.C. Licciardello, and T.V. 

Ramakrishnan, Phys. Rev. Lett. 42, 673 (1979). 
* T.J. Gramila, J.P. Eisenstein, A.H. MacDonald, L.N. Pfeif- 

fer, and K.W. West, Phys. Rev. Lett. 66, 1216 (1991). 
^ R. Pillarisetty, H. Noh, D.C. Tsui, E. Tutuc, E.P. De- 

Poortere, and M. Shayegan, Phys. Rev. Lett. 89 16805 



(2002). 

® Cyclotron resonance measurements on samples A and D, 
give a band mass of 0.17 and O.lSme, respectively. Earlies- 
we had assumed m* — 0.38me. Recent experimental stud- 
ies of the cyclotron resonance in narrow p-type GaAs quan- 
tum wells, similar to our samples, have yielded a band mass 
of 0.19me, which is very close to the values we obtain. See 
W. Pan, K. Lai, S.P. Bayrakci, N.P. Ong, D.C. Tsui, L.N 
Pfeiffer, and K.W. West, Appl. Phys. Lett 83, 3519 (2003). 

^ A. Jauho and H. Smith, Phys. Rev. B 47, 4420 (1993). 



6 



* M. Kellogg, J. P. Eisonstcin, L.N. Pfciffcr, and K.W. West, 
Solid State Conmmn. 123, 515 (2002). 

® E.H. Hwang, S. Das Sarma, V. Braude, and A. Stern, Phys. 

Rev. Lett. 90, 086801 (2003). 
° T.J. Gramila, J.P. Eisenstein, A.H. MacDonald, L.N. Pfeif- 

fer, and K.W. West, Phys. Rev. B 47, 12957 (1993); 
^ H. Noh, S. Zelakiewicz, T.J. Gramila, L.N. Pfeiffer, and 

K.W. West, Phys. Rev. B 59, 13114 (1999). 
^ The barrier in Sample C consists of a 100 A segment of 

AlGaAs, grown in between two 100 A segments of AlAs. 

J.P. Eisenstein, L.N. Pfeiffer, and K.W. West, Appl. Phys. 

Lett. 57, 2324 (1990). 

* For completness, we would like to mention that density 
ratio measurements in the d — 450 A sample exhibit fea- 
tures near matched density, which indicate the presence 
of a small phonon mediated drag contribution. However, 
we believe that the deviation from the dependence in 
this sample arises from the same mechanism that produces 



this deviation in the samples with smaller d. This is since 
we find that the deviations from the dependence at d = 
450 A are unaffected by significantly mismatching the layer 
densities, along with the fact that both the fractional de- 
viation to and the peak position in po/T^ vs T are the 
same as those found at smaller layer spacing. 
S. Das Sarma and E.H. Hwang, Phys. Rev. Lett. 83, 164 
(1999); Phys. Rev. B 61, 7838 (2000). 
G. Zala, B.N. Narozhny, and I.L. Aleiner, Phys. Rev. B 
64, 214204 (2001). 

Y.Y. Proskuryakov, A.K. Savchenko, S.S. Safonov, M. Pep- 
per, M.Y. Simmons, and D.A. Ritchie, Phys. Rev. Lett. 
89, 76406 (2002); H. Noh, M.P. Lilly, D.C. Tsui, J. A. Sim- 
mons, E.H. Hwang, S. Das Sarma, L.N. Pfeiffer, and K.W. 
West, Phys. Rev. B 68, 165308 (2003). 
S. Hani, A. Yacoby, D. Mahalu, and H. Shtrikman, Science 
292, 1354 (2001). 



